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A series of high power, high efficiency A^-band and W-band gyroklystron experiments has been 
conducted recently at the Naval Research Laboratory (NRL). Stagger tuning of the cavities for 
bandwidth enhancement is commonly used in the conventional multicavity klystrons. The desired 
stagger tuning is usually achieved via mechanical tuning of the individual cavities. However, in the 
millimeter wave regime, particularly, in the case of the high average power operation, it is desirable 
to be able to achieve the required stagger tuning by design. The NRL gyroklystron experiments 
explored the tradeoffs between power, bandwidth, efficiency, and gain to study the effects of large 
stagger tuning in millimeter wave without resorting to mechanical tuning of the cavities. Both, 
A a -band and W-band, experiments demonstrated a record power-bandwidth product. The success of 
the experiments was due in large part to a battery of improved large-signal, stability, and cold test 
codes employed in the modeling and design stage. Theoretical models that provide the basis for 
these design simulation tools and the design methodology will be presented. © 7999 American 
Institute of Physics. [S1070-664X(99)92905-1] 
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I. INTRODUCTION 

Gyroklystrons in A a -band and W-band are attractive 
sources for the next generation of millimeter wave radars. 1 
This is because 1 gyroklystrons are capable of operating at 
both high peak power and high average power with high 
efficiency at these frequencies. 2 The instantaneous band¬ 
width of gyroklystrons is determined by the quality factor 
(Q) of the cavities, and especially by the Q of the output 
cavity. This leads to narrow bandwidth, since the Q is typi¬ 
cally on the order of i50-300. However, it is possible to 
increase the bandwidth by employing a multicavity gy¬ 
roklystron circuit, with each cavity tuned to a different fre¬ 
quency (stagger tuning). Stagger tuning of the cavities for 
bandwidth enhancement is commonly used in conventional 
multicavity klystrons. The desired stagger tuning is usually 
achieved via mechanical tuning of the individual cavities. 
However, in the millimeter wave regime, particularly, in the 
case of the high average power operation, it is desirable to be 
able to achieve the required stagger tuning by design. There¬ 
fore, the purpose of our investigations is to develop and vali¬ 
date the simulation models, the design methodology, and the 
manufacturing procedure to achieve this goal for high power 
gyroklystrons. In this paper we describe the simulation mod¬ 
els and the design procedure to address the relevant issues. 

A series of experiments has been conducted recently at 


the NRL to study the effects of large stagger tuning on the 
performance of the gyroamplifiers. These experiments ex¬ 
plored the tradeoffs between power, bandwidth, efficiency, 
and gain due to circuit geometry and operating parameter 
variation. The A^-band experiments, using a three-cavity 3 
circuit with heavy stagger tuning, demonstrated more than 
double the bandwidth of the two-cavity circuit, 4 with only a 
modest decrease in efficiency. The W-band experiments 
achieved a record power-bandwidth product in a four-cavity 
gyroklystron circuit with large stagger tuning. 5,6 For all ex¬ 
periments the results are in excellent agreement with simu¬ 
lations. 

Extensive theoretical work has produced a hierarchy of 
nonlinear models to describe the beam-wave interaction in 
the gyroklystron circuit. The validity of these models and 
their accuracy will be discussed in Sec. II. Ensuring stable 
operation is a central step in the gyroklystron design proce¬ 
dure. Issues relevant to self-oscillations in the drift sections 
of the gyroklystron circuit will be discussed in Sec. HI. The 
design procedure of the input coupler will be described in 
Sec. IV. Finally, in Sec. V the simulation results are com¬ 
pared with the experimental data and the methodology de¬ 
veloped for designing the high average power W-band gy¬ 
roklystron is briefly discussed. 
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II. THEORETICAL MODEL 


The combination of the phase bunching due to the cy¬ 
clotron maser instability in a cavity and the ballistic phase 
bunching in drift tubes makes the gyroklystron capable of 
operating at high efficiency, high power, and high gain. Sev¬ 
eral models describing the gyroklystron interaction have 
been developed over the years, ranging from an analytical 
“point-gap” approach 7 to a multimode, self-consistent ap- 
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proach implemented in a simulation code. 8 A complete de¬ 
scription of the variety of theoretical models, approaches, 
and methods is beyond the scope of this paper. In this section 
we will describe briefly the theoretical models, which have 
been implemented in the large signal codes used at the Naval 
Research Laboratory (NRL) for designing the gyroklystron 
experiments. 

We started from the model described by P. Latham 
et al. 9 The assumptions in that model were dictated by the 
intended gyroklystron application, drivers for the next gen¬ 
eration of linear colliders for particle physics research. The 
collider application required a narrow band, low average 
power gyroklystron capable of producing a very high peak 
power with high efficiency. This was achieved by using a 
moderately relativistic 500 kV and 200 A electron beam. 10 
The NRL gyroklystron development is driven by a radar ap¬ 
plication, which requires operation at high average power 
and wide bandwidth, while maintaining relatively high peak 
power and efficiency. 11 This can be achieved by using a 
weakly relativistic electron beam (50-70 kV and 4-10 A) in 
a stagger-tuned multicavity circuit. 

It turns out that the model in Ref. 9, when augmented 
with additional capabilities required for the design of multi¬ 
cavity gyroklystron with stagger tuning for bandwidth en¬ 
hancement, contains essentially all of the effects needed for 
an adequate description of the beam-wave interaction in a 
gyroklystron. Specifically, a more accurate model of the in¬ 
put coupler and a model for calculating the dielectric loading 
of the intermediate cavities have been added. In addition, it 
turned out that the procedure for designing a gyroklystron 
amplifier used in Ref. 9 is not suitable for effective evalua¬ 
tion of the bandwidth of the circuit. This deficiency was 


removed as well. All of these efforts resulted in a new pro¬ 
cedure for designing wide band gyroklystrons and to an im¬ 
proved version of the large signal code, which we named 
MAGYKL. Here we will describe only the additions to the 
model, keeping the same notation as Ref. 9. 

We start with the evolution of the complex amplitude in 
each cavity in MAGYKL, which is described by a time- 
dependent equation 
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where a is the complex normalized amplitude of the fields; 
A Ct> = (2/{Re{o> c }—is the normalized frequency shift, 
is the load quality factor of the cavity, co is the drive fre¬ 
quency, and (o c is the cold resonant frequency, J is the beam 
current density, W EM is the field stored energy, E c is the cold 
cavity electric field, J is current density, T=2rrl(o is the 
period, and r is the slow time 12 normalized to Qilco. Note 
that in the original formulation (Ref. 9), Eq. (1) does not 
contain the time derivative. We will come to this point later 
when addressing the design procedure. 

The cavities are modeled by a series of straight uniform 
sections with abrupt radial discontinuities at the junctions. 
The fields in each section, expanded as a radial series of TE, 
TM, and TEM modes, are determined through a scattering 
matrix solution. 13 In order to accurately predict the amplifier 
bandwidth of the amplifier, the input cavity model used in 
the past was modified to include a frequency-dependent elec¬ 
tric field amplitude in the drive cavity, as dictated by the 
resonant frequency and total Q of the input cavity. The field 
amplitude in the drive cavity is given by 
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where P m is the input power, which is assumed to be fre¬ 
quency independent, G external * s the external quality factor, 
Glosses is the quality factor which includes ohmic and any 
other losses from the input cavity, such as field leakage 
through the beam tunnel, v 1 and v z are the perpendicular and 
axial electron velocities, £=za)lc is the normalized distance, 
and the angle brackets denote an average over particles. In 
the derivation of Eq. (2), we assumed that the electric field in 
the input cavity is small enough that the beam-wave interac¬ 
tion can be evaluated in the linear regime. 

Before describing the design procedure, a few remarks 
concerning usage of “reduced models” are in order. It was 
originally expected that the model used in Ref. 9 was unnec¬ 
essarily detailed for a weakly relativistic and low current 
beam. Therefore a reduced model 14 based on a single mode 
sinusoidal or Gaussian axial profile for the “cold” field 



could be used. Although reduced models of gyroklystron am¬ 
plifiers are still useful for understanding the basic operation 
of these devices, they are not adequate as a design tool if the 
combination of the highest possible bandwidth and efficiency 
is desired. We illustrate the limitations of the reduced models 
by comparing the results produced by these models with our 
improved large signal code, MAGYKL. We will use the four 
cavity W-band gyroklystron with stagger tuning as an ex¬ 
ample. If the “cold field” in the cavities is represented only 
by one radial mode, TE 01 , and the axial profile is sinusoidal, 
the output power versus input power (the transfer function) 
will be significantly different from the transfer function pre¬ 
dicted by MAGYKL. This is shown in Fig. 1. Similarly, the 
power versus frequency dependence (bandwidth curve) at 
saturation as predicted by magykl is quite different from the 
prediction of a model based on a single mode sinusoidal or 
Gaussian axial profile for the “cold” field. 
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FIG. 1. Efficiency versus input power in the case of four-cavity gy- 
roklystron. The four-cavity structure with the following parameters: Cavity 
1 f c = 93.4064 GHz, <2 ext = 150, £ loss =800; Cavity 2 f = 94.2063 GHz, 
2ioaded= 175; Cavity 3 f c = 93.2836 GHz, Q m<t<x = 175; Cavity 4 f c 
= 93.8874 GHz, Q loaded = 162.43; electron beam parameters are: y 
= 65kV, / = 6A a =1.5, guiding center radius r g = 0.081 cm, zero velocity 
spread. Magnetic field was optimized for maximum efficiency, B 
= 36.8 kG. 


A. Design procedure 

The process of using magykl for designing a gy- 
roklystron is a multi step procedure, which includes the 
choice of an electrodynamical structure, the electron beam 
and drive signal parameters, and the optimization of device 
performance by varying the magnetic field. The determina¬ 
tion of the electrodynamical structure consists of several 
steps. First, the resonance frequency, Q , and length for each 
cavity are chosen. Based on this information, CASCADE 13 is 
used to determine approximate cavity shapes and sizes and 
the drift tube radii that yield the desired frequencies and Q s 
as well as the electric field profile, E c (r,z). The field axial 
profile and the amplitudes and phases of the modes together 
with the “cold” frequencies and Q s are used in the large 
signal code MAGYKL, in which the particle equations and the 
wave amplitude Eq. (1) are solved. At this point the input 
cavity parameters, such as drive frequency, P in , Q ext , and 
<2 loss , the beam parameters, such as current, voltage, 
a~v L lv z , root-mean-square (rms) velocity spread and the 
guiding center radius, and the magnetic field value are cho¬ 
sen. 

The amplitude in the input cavity is determined by Eq. 
(2) via iterations. Usually, three to five iterations are suffi¬ 
cient. The number of the time steps needed for the ampli¬ 
tude, a , for the cavities other than the input to reach a steady- 
state solution is usually on the order of 100. The overall 
efficiency is evaluated using this steady-state amplitude. 
Note that the steady-state solution is consistent with the pre¬ 
scribed total quality factor. 

In the original version of the large signal code, 9 without 
the time derivative term, the values Q and the magnitude and 
phase of the complex amplitude a and its phase and Q were 
used to determine A w from the energy balance equation, as 


in Eq. (7) in Ref. 9. This approach proved to be computa¬ 
tionally effective for the cases studied in Ref. 9. However, it 
turned out to be not practical when stagger tuning was em¬ 
ployed for both efficiency and the bandwidth optimization. 
Adding the time derivative did not change the nature of the 
final solution, since the complex amplitude has to reach a 
steady state in order to satisfy the power balance equation. 
However, the numerical implementation of solving Eq. (2) is 
such that the amplitude and phase of the wave are adjusted 
from one time step to another automatically until the power 
balance equation is satisfied. The computation is repeated for 
different values of the magnetic field in order to optimize the 
performance. Note that the drift length is mainly determined 
by the gain or efficiency optimization. 

The choice of the guiding center radius involves addi¬ 
tional considerations. On the one hand, it is desired to make 
the field of the operating mode evanescently small in the 
drift region, hence a small drift tube radius is advantageous. 
On the other hand, the beam size has to be adequate to en¬ 
sure maximum coupling to the electric field and minimum 
interception by the drift walls. Usually a distance on the 
order of three Larmor radii from the drift wall is considered 
to be safe. Another interesting point worth mentioning con¬ 
cerns the adequacy of the small Larmor radius approxima¬ 
tion that is typically used for describing the beam-wave in¬ 
teraction in the weakly relativistic regime. We find that if the 
small Larmor radius approximation is used, the frequency 
pulling due to the beam loading is about 50 MHz different 
from MAGYKL calculations (no small Larmor radius approxi¬ 
mation is used in Bessel function evaluation in the equations 
used in MAGYKL). By itself it is a tiny error, blit considering 
the fact that the stagger tuning of the cavities is on the order 
of 100-150 MHz, the 50 MHz error is quite significant. 

B. Effects of the dielectric insert 

The total quality factors for the intermediate cavities are 
assumed known in our model. In practice, in order to realize 
the required cavity losses, a ring of a lossy dielectric is in¬ 
serted into the cavities. MAGYKL requires prior calculation of 
the field eigenmode profiles for each cavity, which are nor¬ 
mally obtained from the scattering matrix code CASCADE, 
which does not allow the ceramic insert to be included in the 
computation. Therefore some approximation was necessary 
to model the actual device. Using MAGYKL, a design begins 
with a set of frequency and Q values for each cavity, and a 
known dielectric constant for the ceramic inserts. 

An unpublished code, dielec, developed by P. Latham 
and modified by B. Levush and A. Berman, is used for each 
cavity to find dimensions (radius, length of the dielectric 
ring), dielec is a more advanced scattering matrix code that 
permits (lossy) ceramic rings in the description of the cavity 
geometry. Loss is included through the use of complex val¬ 
ues of the dielectric constant. The cavity frequency and Q 
obtained are passed as inputs to magykl. Note that, when 
DIELEC is used, the number of modes required to converge on 
the eigenfrequency is about 10. However, since cascade 
must be used for the field profiles passed to magykl, ap¬ 
proximate profiles are computed with the assumption that the 
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cavity electromagnetic fields are excluded from the dielectric 
ring as if it were metal. This assumption is valid if the real 
part of the dielectric constant of the ceramic is large. Cavity 
dimensions for this equivalent cavity are kept the same as 
those of the true cavity, except that the radius used in CAS¬ 
CADE is adjusted slightly to match the required resonant fre¬ 
quency. We have recently adapted magykl to take fields 
from dielec. Detailed description of the model will be pre¬ 
sented in other publications. The adapted code, MAGYKL-C, 
follows a similar algorithm, but instead of requiring an ap¬ 
proximate equivalent cavity to be constructed, the field pro¬ 
file is taken directly from DIELEC. The cavity frequency and 
Q are therefore fully consistent with the cavity and dielectric 
dimensions, so that a level of approximation has been re¬ 
moved from the calculation. 

The original code magykl has been used successfully in 
the design of NRL gyroklystron experiments using the cavity 
field approximation described, and therefore for the param¬ 
eter regime of those experiments the approximation of field 
exclusion from the dielectrics seems most likely valid. Here 
we explore this validity by comparing results based on the 
original model with the more complete model, and extend 
the comparison to other parameter regimes of experimental 
interest. We will consider two cases of dielectric properties: 
(A) large dielectric constant, high loss (e= 12.24-/3.67) 
and (B) large dielectric constant, moderate loss (e=10.0 
— /0.88). MAGYKL has been used successfully in the design 
of devices in the case “A” regime, taking the field profile as 
if the ceramic were metal. Case “B” requires a greater 
length of ceramic insert, necessary to achieve the same over¬ 
all absorption to achieve the design Q. The DIELEC field pro¬ 
file (“true” cavity) has a more complex structure, and the 
cascade field profile used by magykl is not as good an 
approximation in this parameter regime. Figure 2 shows the 
profile of the cavity field for the radial position at which the 
electron beam travels, for both the true cavity and the 
equivalent cavity for case B. The effect of the extended di¬ 
electric region is apparent, causing a plateau in the field pro¬ 
file. This deviation from the equivalent cavity profile is suf¬ 
ficient to influence the electron beam. Figure 3 shows the 
transfer function curves calculated by both magykl and 
magykl-C for this case. It can be seen that magykl and 
MAGYKL-C differ distinctly in the peak efficiency. Similarly, 
the efficiency versus frequency dependence as predicted by 
MAGYKL and MAGYKL-C is different. For this design, it 
would be necessary to use the new version of the code to be 
confident of the results. However, in the cases, such as case 
A, the benefits of the computation speed justify the proce¬ 
dure of using th^ “cold” field axial profile in MAGYKL with¬ 
out taking into account the effect of the lossy dielectric. 

III. DRIFT TUBE DESIGN 

Gyroklystron drift spaces ideally provide a field-free en¬ 
vironment in which ballistic bunching of the electron beam 
can occur. As discussed in preceding sections, the drift tube 
radius is made small enough to ensure that the cavity oper¬ 
ating mode is cut off inside the drift tube. However, this does 
not guarantee that the drift tube is free from the fields asso- 
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FIG. 2. The cavities’ parameters are as follows: Cavity 1 f c 
-92.486 GHz, <2ext = 300, <2ioss“800; the drift length between the first and 
second cavity £^1 — 0.64 cm; Cavity 2 length Z, 2 = 0.5cm, f c 
= 93.494 GHz, Qi oa ded“ 170; the drift length between the second and third 
cavity L drift2 =0.82cm; Cavity 3 length L 3 = 0.5cm, f c ~ 92.917 GHz, 
£ loaded = 170; the drift length between the third and forth cavity Ldnfu 
= 0.82 cm; Cavity 4 f c = 93.315 GHz, £i oa ded = 338; the dielectric permittiv¬ 
ity, e = 10.00— iO.88 has a slightly lower real part, permitting greater field 
penetration, and a smaller loss tangent, so that a larger region of dielectric is 
necessary to achieve the same value of design Q. This was achieved with the 
dielectric inserts of width L diel cavl = L diel cav2 =0.18 cm. The solid line in the 
figure is the electric field profile in the second cavity determined by cas¬ 
cade. The dashed line in the figure is the electric field profile in the second 
cavity determined by dielec. 

ciated with self-oscillation. A simple metal drift tube, 
bounded by transitions to the gyroklystron cavities on either 
end, can function as an unintentional resonant cavity to 
modes different from the normal operating mode. This is 
particularly true with gyroklystrons operating in the TE 01 
mode, since the lower order TE n mode can propagate at drift 



P h (W) 


FIG. 3. The transfer curve predicted by magykl (dielectric insert assumed 
not perturbing the fields in the cavity, hence cascade fields are used) and by 
magykl-C (dielec fields are used). The electromagnetic structure parameters 
are the same as in Fig. 2. The beam parameters are: V=55kV, 7=4A, a 
= 1.5, rms (Ai/ X lv L )- 10%. Magnetic field is 35.8 kG and the drive fre¬ 
quency is 93.45 GHz. 
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tube radii and frequencies considerably below the limiting 
values for the TE 0 i mode. These types of instabilities asso¬ 
ciated with inadvertent drift tube resonators, as well as 
backward-wave gyro-instabilities in the drift tubes, become 
increasingly troublesome as one pushes toward higher beam 
power. 

To prevent drift tube oscillations, lossy dielectric mate¬ 
rials (absorbers) are placed inside the drift spaces. 15 The ab¬ 
sorber consists of a continuous liner of lossy dielectric that 
covers most of the drift tube length. The extreme ends of the 
drift tube adjacent to the cavities are metal, to prevent power 
deposition into the dielectric from leakage fields. Often ithe 
dielectric radius is made somewhat larger than the surround¬ 
ing metal, to help prevent charge buildup from stray elec¬ 
trons from the beam. The procedure for analyzing the stabil¬ 
ity of a candidate drift tube consists of a number of steps that 
are applicable to any of the configurations, although the de¬ 
tails will depend on the specific structure under investigation. 
Essentially, one first determines the frequencies where oscil¬ 
lation could occur. This is followed by an evaluation of the 
attenuation and quality factor of the lossy structure at these 

i 

frequencies, in comparison to computed start oscillation val¬ 
ues. Finally, if needed, the physical dimensions of the drift 
tube structure can be adjusted to optimize performance in an 
iterative manner. 

To illustrate the procedure, consider placing a liner with 
a radial thickness of 0.8 mm and a complex dielectric con¬ 
stant of s = 12.24-“ f 3.67 around the 1.75 mm drift tube. Pro¬ 
vided that the length of the dielectric is at least several wave¬ 
lengths long, one can ignore edge effects and model the 
structure as being infinitely long. 

Cyclotron maser instabilities are most likely to occur at a 
frequency where the dispersion curve (a plot of frequency 
versus axial wave number) of the drift structure intersects 
with the dispersion curve for Doppler upshifted cyclotron 
waves supported by the electron beam. When a dielectric 
liner is placed into the drift tube, two effects occur. First, for 
modes that can propagate in the empty guide at the fre¬ 
quency of interest, the presence of losses creates a complex 
propagation factor in the axial direction f3=k z —ik zi , where 
the axial variation of fields is assumed to follow the func¬ 
tional form e~iP z . The term k zi represents an attenuation 
caused by the losses. The second effect of the dielectric is a 
significant change in k z ; modes that were originally heavily 
cut off at the frequency of interest (k zi > 0, k z = 0) can often 
propagate as a result of the addition of dielectric. For the 
lossy dielectric liner to be effective, it is important that these 
newly liberated modes receive sufficient attenuation as well. 
In the same manner, the addition of dielectric also modifies 
the functional form of the dispersion curve and changes the 
instability frequencies. The modified dispersion curves for 
what originally were the TE n , TM U , and.TE 12 modes are 
shown in Fig. 4, and are now indicated by HE n , HM n , and 
HE 12 . This latter notation indicates that the modes are now 
hybrid (both E z and H z components exist). The new HE n 
backward-wave intersection frequency with the beamline is 
now at 51.8 GHz, while the forward-wave interaction with 
this mode is eliminated. However, the HM n and HE 12 curves 
have been pulled down to have intersections at 56.0 GHz for 
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Axial wavenumber k (cm' 1 ) 
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FIG. 4. The dispersion curves are plotted for the TE U , TM n , and TE 12 
modes for a =1.75 mm. The dispersion curve for cyclotron waves on a 
relativistic beam is given* by a)=£l+k z v z , where H is the relativistic cy¬ 
clotron frequency, k z is the axial wave number (wave propagation factor), 
and v z is the axial beam velocity. It is represented by the straight line in the 
figure, for the case of a 65 kV beam with a 36.8 kG magnetic field and a 
= 1.5. The modified dispersion curves for what originally were the TE n , 
TM n , and TE J2 modes are also shown in the figure, and are now indicated 
by HE n , HM n , and HE 12 . This latter notation indicates that the modes are 
now hybrid (both E z and H. components exist). 


the backward wave HM n , 75.2 GHz for the backward-wave 
HE 12 , and 119.2 GHz for the forward-wave HE 12 . Having 
identified the relevant frequencies for beam-wave intersec¬ 
tions (and thus potential instabilities), one can now evaluate 
the stability of the loaded structure. 

To better quantify the effect of dielectric on each of the 
modes, one must realize that is often not the attenuation per 
unit length that is critical for suppressing oscillation, but 
rather the quality factor of the inadvertent cavity that can 
exist within the drift tube. This quality factor will determine 
whether a particular beam power will trigger oscillation 
since stability scales as the product of beam power with Q. If 
the transitions on both ends of the drift tube produced total 
reflection (a worst-case situation), the dielectric-induced 
quality factor of such a cavity would be Q = Qn=k z lk zi . In 
Fig. 5 the quality factors of the modes are plotted versus 
frequency. The potential instability frequencies from Fig. 4 
are indicated with the circles. One can see that quality factors 
at the HE n , HM n , and the lower HE 12 instability frequen¬ 
cies are below 7, which indicates extremely effective load¬ 
ing. Since these values of Qn are massively below the empty 
beam tunnel starting Q for the design current (about 200), 
drift tube stability at these intersection frequencies is guar¬ 
anteed. Because the gyroklystron amplifier will be used in 
radar, it must be unconditionally stable at the operating point 
in the cavities as well as the drift tubes. Thus, a linear analy¬ 
sis was used to determine stability conditions in each part of 
the circuit. 16 More detailed stability analysis is out of the 
scope of this paper and will be published elsewhere in the 
future. 
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FIG. 5. The quality factors, Q = k z lk zi , of the modes are plotted versus 
frequency. The potential instability frequencies from Fig. 4 are indicated 
with the circles. 


IV. DESIGN OF THE INPUT COUPLER 

The input coupler has to satisfy the following require¬ 
ments. First, the bandwidth of the coupler should be enough 
to meet the demands of the gyroklystron bandwidth. Second, 
the operating mode should be TE 0 i and there should be high 
degree of mode purity. Third, the amount of power coupled 
into the TE 01 mode should be maximized. To meet these 
requirements we choose a coaxial slotted structure for the 
input coupler, 17 shown in cross section in Fig. 6. The coupler 
has the advantage of a single guide feed, good cavity mode 
purity, minimal ohmic loading of the cavity walls, and it can 
be designed with a low loaded quality factor. 

The input wave guide excites the TE 41} mode in the 
outer cavity. Power is coupled through four axial slots to the 
inner cavity TE 011 mode. We employed in the design effort 
two general three-dimensional (3D) codes, HFSS 18 and 



FIG. 6. hfss simulation of coaxial drive cavity in cross sectional view 
showing I#! 2 . 
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ARGUS. 19 These codes are the only practical tool to investi¬ 
gate the effects of cylinder/coax centerline offsets, coupling 
slot axial offsets, nonequal slot azimuthal spacing, slot rota¬ 
tions, cavity wall distortions, and other nonideal configura¬ 
tions. These codes are viewed as complementary in their 
characteristics; each with particular advantages that can be 
exploited. The tetrahedral modeling elements employed by 
hfss conform nicely to structures with curved surfaces, al¬ 
though small “facets” must be used for high accuracy. This 
results in a large number of unknowns, leading to large prob¬ 
lems that strain computational resources. Geometric param¬ 
eterization is also difficult due to the adaptive meshing, since 
a geometry change results in a different initial mesh. On the 
other hand, the ability to transmit a wave through the struc¬ 
ture at a particular frequency closely mimics typical testing 
or operation scenarios. ARGUS clearly identifies the cavity 
eigenmodes, in contrast to HFSS, where decomposition of the 
electric field is required to determine mode content. The 
structured grid of ARGUS results in a “stairstep” approxima¬ 
tion of curved surfaces, requiring a dense grid for adequate 
resolution. ARGUS is somewhat faster* since it can be run on 
Cray supercomputers, although the overall solution time is 
comparable to that using HFSS, since multiple runs are re¬ 
quired for Q determination. Large problems can be handled 
with ARGUS, albeit with long solution times. However, pa¬ 
rameterization is relatively easy with a structured grid. Fi¬ 
nally, the real utility of using both codes for coupler model¬ 
ing is confidence in the results. The same answers are 
obtained using two completely different solution approaches. 
Both codes were used with various methods of solution to 
model many issues concerning the input coupler, including 
resonant frequencies, modes, g’s, resistive loading, mode 
purity, parameterization, and manufacturing tolerances. 1 In 
particular, ARGUS and hfss were used to model the Q and 
resonant frequency for the operating mode of the input cou¬ 
pler. Excellent agreement (0.1 %) was achieved between the 
two codes. 

Several techniques were employed with HFSS to deter¬ 
mine the resonant frequency and Q of the TE 01 mode in the 
cylindrical cavity. The cavity resonant frequency and Q were 
determined from the frequency dependence of stored energy. 
The resonant frequency, f r , is the frequency at maximum 
stored energy and Q = f r Kfu~f\)> where f u and are the 
upper and lower frequencies where the stored energy is half 
of the maximum value. Specifically, E*E was calculated 
over the cavity volume, using the “fast-frequency sweep” 
option which uses an asymptotic wave form evaluation to 
calculate a “bandwidth” solution from a full field solution at 
a single frequency. Excellent agreement with ARGUS results 
was obtained with a sufficiently high number of polygon 
segments approximating circular boundaries, ~2°/facet. A 
second method “de-embeds” S n at the drive port to find the 
“detuned short” position. The resonant frequency and Q are 
calculated from the phase of the reflection coefficient. HFSS 
has also been used to model a ‘ ‘transmission’ ’ measurement 
by launching the TE 01 mode in the cylindrical cavity by eva¬ 
nescent coupling through the drift tube and “receiving” at 
the rectangular port. 

ARGUS is used in the eigenvalue solver mode. In prin- 
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ciple, all relevant modes and frequencies can be determined. 
Since the entire problem volume, including conductors, is 
meshed, a sector of the configuration is typically modeled for 
reasonable problem size and solution time. Resonant fre¬ 
quencies and Q s in ARGUS are determined from a variation 
of the “detuned short” method due to Kroll and Yu 20 requir¬ 
ing multiple runs as the length of the input wave guide is 
varied. The method of geometry input in ARGUS is more 
conducive to parametric variation than for HFSS, once the 
basic geometry has been set. HFSS, on the other hand, can be 
more quickly set up when major geometry changes are in¬ 
volved. Both HFSS and ARGUS were used to estimate the 
mode purity of the TE 011 circular cavity mode. Performance 
of the input coupler is found using both HFSS and ARGUS to 
be quite sensitive to dimensions, and the dimensional sensi¬ 
tivity approaches the limits of ordinary manufacturing toler¬ 
ances. Note that a 0.0003 cm change in cylindrical cavity 
radius shifts the resonant frequency by —100 MHz. It is 
important that the slot angular orientation be within about 1° 
of the normal value of 45° with respect to the input wave 
guide. A 1° slot rotation increases/by —20 MHz and re¬ 
duces Q by 3%. A 5° slot rotation increases/by —36 MHz 
and reduces Q by 30%. A frequency and Q within ±100 
MHz and ±50, respectively, of the design values are accept¬ 
able from overall circuit performance considerations. 


Levush et al. 2239 



FIG. 7. Theoretical (solid line) and experimental (black circles) peak power 
and efficiency as a function of drive frequency. Cavity 1 f c = 93.00 GHz, 
G ext =150, <2 loss — 800; Cavity 2/ c = 93.74 GHz, e Ioaded =108; Cavity 3 f c 
= 93.14 GHz, £>i oaded = 107; Cavity 4 / c = 93.49 GHz, 0i oaded =145; Elec¬ 
tron beam parameters are: V=57.72kV, /=4.2A, a=1.5, perpendicular 
velocity spread, Av L /v 1 —9%. Magnetic field was optimized for maximum 
efficiency and bandwidth above 600 MHz, B = 36.0kG. The 3 dB band¬ 
width is 640 MHz. 


V. SUMMARY 

The design of both two and three cavity K a -band gy- 
roklystron circuits have been validated experimentally. 3,4 
The design of the W-band gyroklystron circuits have been 
validated experimentally with record performance in effi¬ 
ciency, peak power, and bandwidth, 5,6 thereby demonstrating 
the validity of the design codes and the design methodology. 
For example, in Fig. 7 (Fig. 8 from Ref. 6) the measured 
efficiency and peak output power versus frequency are plot¬ 
ted together with the predicted values. The theoretical pre¬ 
dictions are in very good agreement with the measured data. 
In this particular example the W-band gyroklystron circuit 
performance was designed to mimic the electrical perfor¬ 
mance of the high average power version of the gy¬ 
roklystron. Radar applications usually require devices having 
high peak and average power capabilities. The procedure of 
the circuit, drift section, and input coupler design described 
above has been used in designing an industrial version of the 
NRL W-band gyroklystron, namely a four-cavity gy¬ 
roklystron capable of operating at 100 kW peak power, with 
efficiency more than 20% and instantaneous bandwidth of 
600 MHz and capable of delivering 10 kW average power. 11 
The 3D codes HFSS and ARGUS have been also used to iden¬ 
tify the geometry of the intermediate cavities that provides 
the assumed eigenfrequencies, quality factors, and mode pu¬ 
rity, while being consistent with thermal management and 
material constraints. Furthermore, magykl was used to 
study variation in circuit performance with changes in circuit 
dimensions, thus establishing manufacturing tolerances. 21 
The thermal management issues for the cavities have been 
addressed by using the 3D thermomechanical code ANSYS. 22 
The performance of the high average power W-band gy¬ 


roklystron will demonstrate the extent of the validity of the 
entire design methodology, which includes design tools for 
both electrical and thermomechanical analysis. 
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